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R.  R.  Putz* 

ABSTRACT 

Tb*  establishment  of  quantitative  relationships  between  recorded 
wave- system  characteristics  and  other  phenomena  requires  numerical  descrip- 
tion of  the  wave  record.  Current  conoepts  applied  to  time  histories  of  wave 
aotivity  at  a point  are  discussed.  Characteristic  statistical  regularities 
found  in  wave  measurements  are  described.  Examples  given  show  the  applica- 
tion of  statistical  techniques  to  the  description  of  wave  systems  in  terms 
of  the  distribution  of  spectral  energy  as  well  as  the  distributions  of 
'•‘individual”  wave  heights  and  periods,  kesults  from  the  prediction,  from 
one  point  to  another,  of  surface  time  histories  Illustrate  the  application 
of  approximate  spectral  information. 

INTRODUCTION 

In  the  study  of  waves  and  their  interaction  with  their  environment, 
there  has  been  a need  for  effective  description  of  observed  wave  systems. 
Convenient  data  (Snodgrass,  1951,  1952)  is  provided  by  the  reoorded  time 
history,  taken  above  some  fixed  point  on  the  bottom,  of  the  fluctuating 
pressure  below  the  surfaoe  of  of  the  water  level  at  the  surface.  It  will 
be  understood  that  the  length  of  a wave  record  selected  for  analysis  is 
long  relative  to  the  periods  of  the  fluctuations  of  interest,  but  short 
compared  to  the  variations  in  meteorological  conditions. 

The  questions  we  are  concerned  with  are*  ;A\at  information  is  con- 
tained in  a wave  record?  And  how  may  this  information  be  conveniently  ex- 
traoted?  We  shall  first  desoribe  the  statistical  regularities  suggest- 
ing a particular  mathematical  model  (Tukey,  1950;  Pierson  1952a)  found 
useful  for  interpreting  wave-record  data.  The  two  complementary  aspeote 
of  this  model  each  involve  a distribution  - one  a statistical  distribution 
of  wave-rooord  ordinates,  the  other  a speotral  distribution  of  energy. 
Theoretioal  relations  connecting  various  model  parameters  will  be  illus- 
trated with  aotual  data.  Although  the  experimental  results  to  be  pre- 
sented have  been  obtained  with  pressure  waves  in  the  ooean  for  the  most 
part,  many  of  the  methods  and  results  are  applicable  to  other  kinds  of 
data. 


In  the  lower  half  of  Fig,  1 is  shown  a short  segment  of  a typical 
wave  reoord,  giving  tho  pressure  fluctuation  64  feet  below  the  surfaoe 
off  tho  California  ooast.  ^ie  total  time  interval  shown  represents  about 
four  and  one-half  minutes,  during  which  from  twenty  to  twenty-five  waves 
pass  the  reoording  point.  It  will  be  seen  that  on  the  time-history  curve 
the  apparent  slope  and  curvature,  as  well  as  the  ordinate,  vary  irregular- 
ly from  one  instant  to  the  next  and  bear  little  obvious  relation  to  one 
another.  Tho  most  noticeable  feature  may  be  the  tendency  for  curvature 
and  ordinate,  measured  from  its  average  .level,  to  have  opposite  signs, 
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The  upper  half  of  »‘ig4  1 shews  a sample  segment  generated  by  the  theore- 
tics.! mod«l  in  which  suitable  parameters  wore  chosen  to  yield  an  artifi- 
cial wave  record  comparable  to  the  observed  one. 

DISTHIHJ TIONS  OP  ORDINATES  AND  DERIVED  QUANTITIES 

Suppose  a horizontal  line  is  drawn  cutting  the  wave-reoord  ourvo 
at  an  arbitrary  level.  The  fraofcion  of  the  time  the  curve  spends  b6lovr 
such  a line  is  a function  (Birkhoff  and  Kotig,  1955;  Puts,  1953b)  in- 
creasing with  the  height  of  the  line  knowa  as  the  distribution  function 
for  the  our  vs.  A plot  of  this  function  for  a typioal  twenty-minute  wave 
reoord  is  shewn,  in  Pig.  2.  Here  the  number  on  the  vertical  scale  repre- 
sents the  arbitrary  chart  level,  the  horizontal  scale,  the  percent  of  time 
that  the  ovavo  lies  bdlow  that  level.  The  notion  of  the  probability 
of  finding  the  curve  below  a given  level  may  bo  introduced  if  one  thinka 
of  choosing  at  random  an  instant  of  time  on  the  ohart.  The  location  of 
the  plotted  points  on  an  approximate  straight  line  is  characteristic 
(PutJB,  1953b)  and  results  here  from  the  choice  of  the  horizontal  scale 
wrhioh  is  tho  familiar  Gaussian- distribution  or  normal-probability  scale. 

The  interesting  thing  about  wave  records  is  the  apparent  ability  of 
the  normal  probability  scale  to  rectify  very  nearly  not  only  the  distri- 
bution curves  for  the  ordinate  ana  tne  first  derivative,  which  have  been 
experimentally  checked (Rudnick,  1951;  Putz,  1953b),  but  also,  according 
tc  the  theory,  the  distribution  curves  for  derivatives  of  all  orders.  In 
eaoh  oase  the  mean  value  below  v/hioh  the  curve  or  derivative  spends  ju3t 
fifty  peroent  of  its  time  will  be  zero.  For  the  ordinates,  this  will  be 
true  because  of  our  convention  of  measuring  them  from  their  mean.  Since 
in  this  way  one  point  is  fixed  on  each  curve,  the  straight-line  plots  for 
the  various  derivatives  of  the  wave-record  will  differ  among  themselv3S 
only  in  their  slopes.  A convenient  measure  of  the  slope  is  the  difference 
between  the  height  of  the  line  at  the  64.1  percent  level  and  its  height 
at  the  60  percent  level.  This  diff^nc*,  called  the  standard  deviation, 
or  root-mean-square  (r.m.  s.)  ordinate,  when  measured  for  tho  k^b  deriva- 
tive, is  denoted  by  . 

The  ordinate  distribution  concept  may  be  extended  by  considering 
iL^re  than  one  ordinate  at  a time.  Given  n arbitrarily-selected  chart 
levels,  the  n-dimensional  generalization  (Cramer,  1946)  of  the  Gaussian 
distribution  is  then  applicable  to  the  fraction  of  the  time  that,  simul- 
taneously, n ordinates,  chosen  at  different-  time  instants  on  the  wave  reoord 
will  lie  below  their  corresponding  chart  levels.  Tho  assumption  that  for 
each  finite  n,  the  selection  of  n ordinates  actually  resuxts  in  a multi- 
dimensional Gaussian <ii 5 "ributioo  Ns. kncWm  as  tho  co-normad  hypothesis-. Suoh 
tribe  idea  is  characterized  by  a set  of  parameters  known  as  the  covariances 
(Cromer,  1946).  These  parameters  may  be  i easily  interpreted  in  terms  of 
the  degree  to  which  th  . various  ordinncoe  do  tormina , or  are  determined  by, 
each  other.  Computation  (Rice,  1914-5 ' tells  us  that  the  ojvarianoe  of 
any  two  quantities  so  distributed  is  proportional  to  coj(  T7  • p ), 
where  P is  the  probability  that  the  two  quantities  have  opposite  algebraic 
signs. 


CUMULATIVE  DISTRIBUTION  FUNCTION  Of  THE  ORDINATES 


A very  convenient  property  qf  these  multiple  ordinate  diatritutions, 
find  one  which  is  physioally  reasonable  except  for  relatively  long  reoorde, 
ia  the  stationary  property  (Lee,  1949;  Puts  1953a),  namely  that  the  dis^ 
tributica  of  any  set  of  ordinates  depends  only  upon  the  differences  in  time 
between  the  corresponding  absoissas.  We  tnen  have  distributions  with  equal 
oovarianoes  whenever  the  corresponding  two  ordinates  are  separated  by  the 
■nma  time  1*2,  T - 1:9:,  the  oovnrienoe  v between  eny  two  ordinates  *e— 
looted  at  ‘times  T and  t -t  r depends  only  upon  r . The  oovarianoe  ia 
thus  a function  “of  the  lag  r , the  funotion  y ( r ) being  known  as  the  co- 
variance function*  (Mann,  1955),  ana,  when  divided  by  its  maximum  value 
y ( o),  as  the  correlation  funotion,  p (r  ),  the  latter  being  the  familiar 
product-moment  correlation  ooeffioient.  Once  assumed  for  the  ordinates, 
the  stationary  and  multinormal  properties  follow  for  derivatives  of  all 
orders. 


WAYS  HP  I GETS 


Interpretations  for  many  of  the  parameters  of  these  underlying 
probability  distributions  may  be  found  in  the  distributions  of  certain 
quantities  which  nay  be  graphically  moasured  on  the  wave  record. 

If  Njj-  denotes  the  number  of  times  the  k^*1  derivative  passes  through 
the  zero  level,  and  if  I is  the  total  length  of  the  wave  reoord  in  seconds, 
then  the  average  number  of  zero-level  crossings  (up  or  down)  per  seoouu  will 
be  Njc/T,  A natural  definition  of  the  mean  (undireotional)  zero-orossing 
frequenoy  is  then  fk  • Nk/(2T)  in  oycles  per  second,  or  o>  k * ZttE^/(ZT)  m 

irNk/T  in  radians  per  second,  Theory  prediots  that  oj  k = er  k ^ 
which  serves  to  relate  the  r,m.  8.  values  of  the  ordinate  and  of  the  first 
and  second  derivatives  to  the  observed  mean  zero-crossing  frequencies  of 
the  ordinate  and  the  first  derivative,  when  k io  taken  to  be  zero  and  one* 

Further  relations  appear  if  we  consider  simultaneous  values  of  or- 
dinate and  second  derivative,  Ihe  value  of  the  ooeffioient  of  correlation 
between  these  two  quantities  is  given  by  p Q s oos  (tt-Po),  where  F0  is 
the  fraction  of  the  time  that  they  have  opposite  signs.  Theory  (Rioe,  1944-8) 
prediots  further  that  f 0 5 -(  w o/  w 1)  = -(K0Al).  A satisfactory 
experimental  oheck  of  the  corresponding  relation  oos  ( 7T-P0)  B -(No/Hq) 
has  beenobtained  by  locating  points  of  infleotion  on  a few  twenty-minute 
records. 


We  may  now  oonsider  th6  distribution  of  wave  heights.  For  this 
purpose  ?re  shall  take  as  fundamental  the  heights  of  the  maximum,  points 
on  the  wave  record  curve,  is«;,  the  points  at  which  the  second  deriva- 
tive is  negative  and  the  first  derivative  is  zero.  The  theory  applies 
also  to  minimum  points  if  the  ymva-record  curve  is  reflected  about  its 
mean  level.  If  these  peak  Hei  measured  as  directed  distances  from 

the  mean  level,  are  averaged,  the  resulting  quantity  p ^ is  just  one-half 

p tv  (Puts,  ib5Ea). 


L’frG 

XX 

P H s(2f)£(-  P o)  & o » illustrated  for  twenty- 


the  ft  V A yA  ij  Q*)l f*  tirpnA 

retioal  relation  is 
minute  wave  reoords  in  Fig,  3,  where  p u . measured  in  chart  divisions, 
is  plotted  against  (NoAl)  a oi  also  in  chart  divisions,  The  straight 
line  shown  has  the  slope  (2w)o  given,  by  the  theory.  It  is  seen  that  the 
average  of  the  so-called  individual  wave  heights  (Putz,  1953b)  can  provide 
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a good  estimate  of  the  r.m.s.  ordinate  rhen  combined  with  a oo..at  of 
the  number  of  times  the  ordinate  and  the  slope  pass  through  zero. 

The  entire  distribution  of  peak  (or  trough)  heights  (Rice  1944-6) 
given  by  the  theory,  depends  only  upon  cr0  and  pQt  and  ic  shown  in  Fig. 

4.  The  vertical  scale  corresponds  to  the  chart  level  measured  from  the 
mean  in  <T  Q-  units  equal  to  the  r.m.s.  ordinate.  The  percent  of  peak 
heights  not  exceeding  this  vertical  level  appears  on  the  horizontal,  the 
Beale  in  this  case  being  chosen  to  correspond  to  the  so-called  Rayleigh 
distribution  (Longuet-Higgiaa  1952;  Lawson  and  Uhlenbeok;  1950,  Knudtzon, 
1949)  which  the  peak  heights  follow  more  and  more  closely  as  the  param- 
eter pQ  tends  to  the  value  minus  one.  It  may  be  observed  that  the  pro- 
bability of  a low  peak  situated  below  the  mean  level  is  just  jg  (i+/>o). 
which  tends  to  zero  as  p 0 tends  to  minus  one,  corresponding  to  a rela- 
tively narrow-band  spectrum. 

Pig.  6.  shows  a typical  observed  distribution  of  peak  heights  for 
& twenty-minute  wave  record,  The  straight  lin9  represents  the  asymptotic 
Rayleigh  distribution,  while  the  curve  represent:  the  distribution  of  peak 
haightG  actually  predicted  from  measured  estimates  ofm0  and  p 0,  ae 
well  as  of  the  mean  ordinate  fJ-0,  The  vertical  level  is  indicated  in 
chart  divisions  on  the  left  and  in  a0“units  on  the  right. 

ENERGY  DISTRIBUTION 

The  statistical  analysis  of  the  periods  between  the  zero  cros- 
sings, or  between  the  peaks  on  a wave  record  is  olosely  related  to  the 
distribution  of  energy  over  the  frequency  or  period  spectrum.  This  as- 
signment of  to  a continuum  of  superimposed  elementary  waves  ia 

specified  by  the  so-oallod  spect-al  distribution  funotion  or  speotrum. 

The  speotrum  provides  a means  of  deriving  and  interpreting  the  parameters 
of  the  distributions  of  wave-record  ordinates  and  derivatives  entering 
into  the  distribution  of  wave  heights. 

The  distribution  of  energy  is  of  fundamental  importance,  of  course, 
for  the  study  of  the  generation,  propagation  and  the  effects  of  waves, 
■While  information  regarding  the  frequenoy  of  occurrence  of  vreive  heights 
or  periods  of  given  magnitudes  appears  occasionally  to  b'  _i  ootly  useful 
for  design  purposes,  a knowledge  of  the  spectrum  oan  be  used  whenever  the 
frequency-response  funotion  oorresponding  to  a general  filtering  action  ia 
’ known. 


The  usual  spectrum  analysis  of  a given  wave-record  sample  results 
iu  the  assignment  of  an  ampxituus  and,  in  addition,  a phase  angle  to  each 
of  a certain  set  of  frequencies.  Ho.verer,  the  direct  dependence  of  suoh 
phase  information  on  the  time  origin  to  which  the  speotru"1  analysis  ia 
attached,  makes  it  unsuitable  as  a possible  parameter  i.  v n stationary 
random  process.  Thus  it  is  natural  to  take  as  fundamental  the  remaining 
napoo  t of  the  spectrum  --  i.e.,  the  amplitude,  or  equivalently,  the  square 
of  the  amplitude,  the  so-called  power  spectrum  (Rice,  1944-5). 


There  is  a funotion  obtainable  from  a sample  wave  record  known  ae 
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FIGURE  5 


PERCENT  OF  WAVE  RECORD  PE  AKS  NOT  EXCEEDING  »i  (or  M) 


the  sample  covarianoe  function  whioh  has  for  its  speotrum  ,iu3t  the 
powor  spectrum  of  the  wa?e  record  itself*  In  many  oase3  the  power  spec- 
trum of  a gi^en  sample  wave  record  is  most  naturally  obtained  by  first 
oomputing  the  covariance  function  for  the  wave  record  v/hicn  was  observed. 
However,  any  problem  whose  solution  involves  the  prediction  of  tho  pos- 
sible sample  wave  reoords  which  might  have  been  observed  before  or  after 
the  given  one  will  naturally  lead  to  the  association  of  a probabilistic 
model  with  each  sample ,wave  record.  This  model  may  be  thought  of  as 
comprising  a whole  sot  of  possible,  infinitely  long,  wave  reoords,  eaoh 
having  the  same  power  spectrum  or  covarianoe  funotion.  The  identifica- 
tion of  this  common  oovariance  function  with  that  desoribed  earlier  in 
connection  with  stationary  multinomial  ordinete  distributions  is  the 
step  which  supplies  the  relation  between  those  probability  distributions 
and  the  spectral-energy  distribution.  The  specification  of  the  qo- 
varianoe  function,  together  with  the  stationary  multinormal  assumption, 
defines  a sc— i&lled  stationary  Gaussian  random  (or  stochastic)  funotion 
fc>r  process).  (Mann,  1953 j Lee,  I960).. 

The  ability  of  suoh  a theoretical  model  to  describe  observed  wave 
reoords  is  illustrated  by  comparing  a sample  taken  from  suoh  a random 
process  with  a typioal  wave  record.  The  upper  curve  in.  Pig.  1 was  gener- 
ated after  first  specifying  three  fixed  numDers  defining  the  spectrum  and 
then  consulting  a table  of  random  numbers  (Wold,  1946). 

If  the  mathematical  moments  of  the  spectrum  exist  (puts.  1953a)  , 
oartain  of  these  can  be  shown  to  be  equivalent  to  the  r,m. a.  values  of 
the  derivatives  of  the  original  wave  record,  while  certain  ratios  of  these 
moments  are  essentially  the  coefficients  of  oorrelation  between  the  deriva- 
tives. Hius  the  r.m. s.  ordinate  crQ  is  the  zero-order  moment  of  the  spoo- 
trum  or  the  total  speotral  energy,  while  p0  is  related  to  the  second- 
order  speotral  moment,  or  more  precisely,  to  a oertain  measure  of  the 
relative  spectral  bandwidth  when  the  latter  is  6mall. 

Ytf.VE  PERIODS 

It  has  been  seen  that  the  of  time  at  which  the  wave  record 

passes  through  Z0**o  or  through  a peak  enter, by  way  of  the  relative  speotral 
bandwith,lnto  a relation  between  the  wave  heights  and  the  r.m.s.  ordinate. 
For  the  analysis  of  the  sps otrum,  the  differences  between  two  successive 
sero-orossinge  may  be  considered  the  fundamental  individual  wave  periods, 
A,ooording  to  the  theory,  for  a long  record  the  list  of  thes6  individual 
wave  periods  in  the  order  of  their  ocourrenoe  is  completely  equivalent  to 
the  spectrum.  Peak- tc- trough  or  t rough- to-peak  time  intervals,  corres- 
ponding to  the  periods  between  successive  zero-crossings  of  the  wave-record 
derivative,  while  more  sensitive  to  the  influence  of  noise,  are  capable  in 
principle  cf  yielding  the  same  information.  On  the  other  hand,  peak-to- 
peak  or  trough-to-trough  time  intervals  (Snodgrass,  1952),  corresponding  to 
the  periods  between  alternate  zero-crossings  of  the  derivative,  would  be  ex- 
pected to  oontain  less  information. 

The  significance  of  the  zero-crossings  for  our  model  lies  in  the  fact 
that  the  abscissas  of  these  points  soparate  the  intervals  during  whioh  the 
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VIS.V9  record  has  constant  sign*  Using  those  abscissas,  we  may  estimate 
the  correlation  function  for  a given  wave  reoord  by  first  determining  the 
fraotion  ft  (r)'of  the  time  that  the  original  wave  reoord  e.nd  the  wave 
reoord  advanoed  by  a time  interval  j have  opposite  sign.  It  is  frees 
the  correlation  function  P(  r ) » cor  ~[7T  P{  r )]  that  the  spectral  in- 
foraabion  is  obtained* 

Pig.  6 shows  a comparison  of  the  estimate  of  p{  r ) obtained  for  a 
twenty-minute  reoord  by  measuring  zsro-orossing  abscissas  and  the  ordinary 
©eiiioate  (Tukey,  1C4.9?  PierGan^e  JUrfe5')  obtained  by  measuring  ordinate? 
evei  ’ second.  This  agreement  between  correlation  funotions  has  been  found 
sufficiently  olose  to  result  in  a relatively  small  difference  between  the 
corresponding  estimates  of  the  wave-record  power  speotrum  from  the  two 
methods.  It  may  be  observed  that  in  addition  to  measuring  1200  ordinates 
on  the  wave  reooru,  nearly  1200  multiplications  and  1200  additions  are  re- 
quired to  obtain  eaoh  point  of  the  solid-line  ourve,  while  for  e&oh  point  of 
the  dotted  line  ourve  it  is  sufficient  to  measure  about  200  absoiscaa  on 
the  wave  reoord  and  perform  about  400  additions  and  200  subtraotione.  Tfhilc 
sithbr  method  leads  to  a substantial  saving  in  labor  over  q straightforward 
numerioal  Fourier  analysis  of  the  original  wave  reoord  (James,  Niohola 
and  Phillips,  1947);  the  zero-oroseing  method  for  the  correlation  function, 
requiring  essentially  no  multiplication,  may  be  carried  out  with  only  an  ai- 
ding machine. 

The  degree  of  reliability  of  the  zero-crossings  exhibited  by  this 
one  experimental  check  on  a ‘ twenty-minute  wave  reoord  may  be  of  some  theo- 
retical interest.  The  close  correspondence,  particularly  for  small  time 
lags,  seen  in  Fig,  6 is  additional  evidence  for  the  general  mathematical 
model  tvs  have  described,  and  for  the  multinormal  property  in  particular. 

It  constitutes  evidence  also  for  the  sampling  adequacy  of  a eve  reoord 
twenty  minutes  in  length.  The  indications  are  that  the  so-oalled  individu- 
al wave  periods,  suitably  defined,  can  furnish  reliable  information  about 
the  Fourier  speotrum  of  the  wave  record. 

This  kind  of  approximate  speotral  information  has  been  used  to  pre- 
dict, by  a least-squ\res  method  (Putz,  1952b),  the  behavior  of  ocean  waves 
at  one  time  and  plaoe  from  their  behavior  at  a time  and  plaoe  nearby. 

The  resulting  prediction  is  shown  in  Fig.  7,  where  the  upper  plot  repre- 
sents the  time  history  of  the  observed  surfaoe  elevation  at  a point  1570 
feet  seaward  of  the  point  at  whioh  prediction  was  to  be  made.  The  lower 
pare  of  the  figure  shows  the  observed  and  the  predicted  wave  record  at 
•the  shoreward  point  fifty  seconds  later.  The  fact  that  the  accuraoy  of 
prediction  represented  here  oompares  favorably  with  that  based  upon  a 
more  exact  spectrum  analysis  is  a further  indication  of  the  reliability  of 
the  information  provided  by  the  zero-crossings. 

SUMMARY 

There  exists  An  approximate  mathematical  model  for  wave  records 
which  may  be  profitably  exploited  until  refinements  in  it  are  found  ne- 
cessary* This  model  specifies  simple  verifiable  relations  which  may  be 
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espeoted  to  hold  between  appropriate  quantities  easily  measured  on  long 
wave  records,  Ihe  parameters  coourring  in  the  model  may  be  estimated  from 
individual  wave-height  and  wave-period  distributions,  A few  such  para- 
meters, easily  interpreted  in  terms  of  the  spectrum,  are  available  for  the 
description  cf  a wave  reoord. 
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